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Abstract 

Exchange interactions in spin systems can give rise to quantum entan- 
glement in the ground and thermal states of the systems. In this paper, we 
consider a spin tetramer, with spins of magnitude i, in which the spins in- 
teract via nearest-neighbour, diagonal and four-spin interactions of strength 
Ji, J2 and K respectively. The ground and thermal state entanglement 
properties of the tetramer are calculated analytically in the various limiting 
cases. Both bipartite and multipartite entanglements are considered and a 
signature of quantum phase transition (QPT), in terms of the entanglement 
ratio, is identified. The first order QPT is accompanied by discontinuities 
in the nearest-neighbour and diagonal concurrences. The magnetic proper- 
ties of a S = I AFM polyoxovanadate compound, V12, are well explained 
by tetramers, with J 2 = 0, K = 0, in which the spins interact via the 
isotropic Heisenberg exchange interaction Hamiltonian. Treating the mag- 
netic susceptibility \ as an entanglement witness (EW), an estimate of the 
lower bound of the critical entanglement temperature, T c , above which the 
entanglement between two individual spins disappears in the experimental 
compound, is determined. Two other cases considered include the sym- 
metric tetramer, i.e. tetrahedron ( J\ = J2, K=0 ) and the symmetric 
trimer. In both the cases, there is no entanglement between a pair of spins 
in the thermal state but multipartite entanglement is present. A second 
EW based on energy provides an estimate of the entanglement temperature, 
Te, below which the thermal state is definitely entangled. This EW detects 
bipartite entanglement in the case of the tetramer describing a square of 
spins ( the case of V12 ) and multipartite entanglement in the cases of the 
tetrahedron and the symmetric trimer. 
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I. Introduction 



Entanglement is a fundamental property of quantum mechanical systems and 
gives rise to an excess of correlations in a system over and above those expected 
from classical considerations [T]. A pure state is said to be entangled if it does 
not factorize , i.e., cannot be written as a product of individual wave functions. 
A well -known example of an entangled state is the singlet state of two spin- 
\ particles, 775 (|ti) ~~ lit))? which cannot be written as a product of the spin 
states of individual spins. Measurement on one component of an entangled pair 
fixes the state of the other implying non-local correlations. In the case of a 
mixed state, entanglement occurs if the density matrix is not a convex sum of 
product states. The importance of entanglement derives from its essential role 
in applications related to quantum information and communication. Candidate 
systems for implementing the application protocols include spin systems in which 
exchange interactions give rise to entanglement EH IH E| • 

Entanglement is a resource which can be created, manipulated and destroyed. 
It can be of different types, e.g., bipartite, multipartite, localizable pj, zero- 
temperature, finite-temperature etc. for which appropriate quantification mea- 
sures are available. Bipartite (multipartite) entanglement involves two (more 
than two) subsystems. The entanglement between a pair of spins belonging to 
a chain of interacting spins provides an example of bipartite entanglement. Bi- 
partite and to a lesser extent multipartite entanglement properties of a variety 
of spin models have been studied so far at both zero and finite temperatures and 
including an external magnetic field [Zl IHl O HOI CH [T21 HBl H31 HSl [TEl - These 
studies show that the amount of entanglement can be changed by changing the 
temperature T and/or the external magnetic field. Since entanglement involves 
non-local correlations of purely quantum origin, an issue of considerable interest 
is whether entanglement develops special features in the vicinity of a quantum 
phase transition (QPT). A QPT occurs at T = and is brought about by tuning 
some system parameter, say, the exchange interaction strength or an external 
variable like the magnetic field to a critical value [TJj. In a QPT, the ground 
state of the system undergoes qualitative changes which in turn affects the en- 
tanglement properties in the ground state. Some recent studies have explored 
the relation between entanglement and QPT in a variety of spin models and the 
main conclusion is that certain entanglement-related quantities exhibit features 
like scaling and singularity in the vicinity of a quantum critical point (QCP) 
H3 DH CHI O EDI EH E2|. In the case of first-order QPTs, the ground state 
concurrences may change discontinuously at the transition point O EH ESI 26j . 
The influence of quantum criticality extends also to finite temperatures so that 
measurements of appropriate observables provide signatures of QPT. At finite T, 
the system in thermal equilibrium is described by the density operator, p (T) = 
j^exp (— ^f), where H is the Hamiltonian, Z the partition function and fee the 
Boltzmann constant. A thermal state remains entangled upto a critical tempera- 
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ture T c beyond which the state becomes separable, i.e., the entanglement falls to 
zero. Experimental demonstrations of entanglement are mostly confined to the 
microworld, i.e., to systems consisting of a few photons, atoms or ions. There 
is now experimental evidence that entanglement can also affect the macroscopic 
properties of solids. This has been shown in the insulating magnetic compound 
LiHoxYi-xFi the specific heat and the susceptibility data of which can only be 
explained if quantum entanglement of the relevant states is explicitly taken into 
account [23 [2H]. Measures of thermal entanglement based on the thermal den- 
sity matrix require a knowledge of both the eigenvalues and the eigenvectors of 
H . On the other hand, there are suggestions that macroscopic thermodynamic 
observables can serve as entanglement witnesses so that a measurement of these 
quantities can provide the evidence for entanglement [HI EHl E3 023- EH- An 
entanglement witness (EW) is an observable the expectation value of which is 
positive in unentangled, i.e., separable states and negative in entangled states 
[H21 EH E3] • The thermodynamic observables which have been proposed as EWs 
include internal energy and magnetization and magnetic susceptibility |281 13 1|. 
The latter has been used as an EW in the spin-^ alternating bond antiferro- 
magnet Cu{NO^)22..hD20{CN). The compound can be considered as a chain 
of uncoupled spin dimers since the ratio of the inter-dimer to the intra-dimer 
exchange interaction strengths is approximately 0.24, i.e., low. For separable 
(unentangled) states, the magnetic susceptibility obeys the inequality 

where g is the Lande splitting factor, fig the Bohr magneton and N the number 
of spins in the system. Entanglement is present in the system if the inequality in 
(1) is violated. The intersection point of the curve representing the EW (equality 
in (1) ) and the experimental % versus T curve defines the critical temperature 
T c below which entanglement is present in the system. The experimental esti- 
mate of T c ~ 5K is in good agreement with the theoretical value of the critical 
temperature at which the pairwise thermal entanglement (entanglement between 
two spins ), as measured by the concurrence, falls to zero. 

Determination of the entanglement properties of an interacting spin system is 
a theoretical challenge as the eigenstates and eigenvalues are not known exactly 
when the number of spins is large. Most of the calculations are confined to systems 
containing a few spins so that exact diagonalization is possible. Studies on finite 
quantum spin systems acquire significant relevance in the context of molecular 
or nanomagnets. In such magnetic systems, the dominant exchange interactions 
are often confined to small spin clusters. The inter-cluster exchange interac- 
tions are much weaker in comparison so that the compounds can be assumed to 
consist of independent spin clusters. A recent study provides a number of exam- 
ples of molecular magnets the thermodynamic and neutron scattering properties 
of which can be well described by small spin clusters like dimers, trimers and 
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tetramers [HSj- As in Ref. [2Hj, one can study the entanglement properties of 
the molecular magnets by treating the susceptibility x as an EW. The earlier 
work dealt with spin dimers for which only pairwise entanglement is possible. 
In this paper, we consider clusters of three (trimer) and four (tetramer) spins 
in which pairwise entanglement between individual spins does not exhaust the 
total entanglement. The tetramer Hamiltonian contains both bilinear and four- 
spin interactions. The ground state and thermal entanglement properties of the 
tetramers are determined analytically. The influence of multispin interactions 
on entanglement is further determined. The system exhibits QPTs at special 
values of the exchange interaction strengths. A signature of the QPT via the 
so-called entanglement ratio is identified. A distinct signature of first order QPT 
is provided by jumps in the amounts of entanglement associated with n.n. and 
diagonal spin pairs. The magnetic properties of the polyoxovanadate compound, 
(NHEt) 3 [V 8 IV VlAs 8 O 40 {H 2 0)\ .H 2 (designated as V12 ) are well explained 
by spin— | AFM tetramers, with only nearest-neighbour (n.n.) interactions, and 
described by the isotropic Heisenberg exchange interaction Hamiltonian [35, 36J. 
The experimental data on the magnetic susceptibility of this compound are avail- 
able. Treating \ as an EW, the critical entanglement temperature, T c , below 
which entanglement is certainly present in the system, is determined. The cases 
of the S = 2 AFM symmetric trimer and tetrahedron are also considered. 

Dowling et al. have introduced the concept of the entanglement gap, de- 
fined to be the difference in the energies of the minimum energy, E se p, that a 
separable state may attain and the ground state energy E . If the energy of the 
system falls within the entanglement gap, the state of the system is entangled. 
The entanglement gap temperature, Te, is defined to be the temperature at which 
the thermal energy U(Te) = E sep , the minimum separable energy. Below Te, the 
thermal state of the system is bound to be entangled. We obtain an estimate 
of T E in the cases of a single square of spins (the case of V12), a tetrahedron 
and a symmetric trimer. In the last two cases the critical entanglement temper- 
ature T c , determined by using \ as an EW, is identical to the entanglement gap 
temperature Te- 



II. Entanglement properties of S = \ AFM tetramer 

We consider a tetramer of spins of magnitude \ (Fig. 1) described by the AFM 
Heisenberg exchange interaction Hamiltonian 

H = J\ [S\.S<2 + 5*2. S3 + 5*3. 5*4 + S4.1S1) + J2\S\.S 3 + S 2 .S&) , . 

+K 1 (S 1 .S 2 )(S 3 .S A ) + K 1 (S 2 .S 3 )(S 1 .S 4 ) +K 2 (S 1 .S 3 )(S 2 .S i ) [Z) 

where Si is the spin operator at the ith site of the square plaquette, J\ is the 
strength of the n.n. exchange interaction, J 2 that of the diagonal exchange in- 
teraction and Ki, K 2 are the strengths of the four-spin exchange interactions. 
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Figure 1: A tetramer of spins of magnitude |. J\ and J 2 denote the strengths of 
the n.n. and diagonal exchange interactions. The four-spin interactions are not 
shown. 

The entanglement properties of the four-spin cluster have earlier been studied 
analytically only for J\ 7^ and numerically for both J\ 7^ 0, J2 7^ [38J. We 
now determine the entanglement properties analytically for the general case in 
Eq. (2). The ^-component of the total spin, 5*°*, is a conserved quantity so that 
the eigenvalue problem can be solved in the separate subspaces corresponding 
to the different values of S 1 " 1 . The results are displayed in the following (Ei , 
i = 1, 16, is the energy eigenvalue): 



'tot 
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^ = ITTTT) 
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^4 = -71 (ITTTI) + ITTIT) + ITITT) + IITTT)) 
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^5 = -71 (ITTTI) + ITITT) - ITTIT) - IITTT)) 
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^ = ^(irur>-um» 

£ 7 =-(f+f+^) (9) 

^ = ^(irm>-imr» (m , 

TP — ( T i J2 3Xi , K 2 \ l iU i 



^ = ^ ami) + mm + mm + mm + mm + mm) 



^9=U+f + f + § 



(11) 



^10 = 7i (ITTU) + |UTT> - ITIIT) - IITTD) 

P _ / 3J 2 , 3Xi , 9fC 2 A 
^10 - ^ 2~ + T" + 



^11 = 7t5 ( 2 ITITI) + 2 IITIT) - ITTU) - ITIIT) - IIITT) - IITTD) 



E u = (-2 Jx + f + Zfi + f 



(13) 



= -1 



^i2 = ^(UUT>-UTU» 

TP — f J 2 , Ki ,3K 2 \ 

iks = ^(IUU>-|TW» 

TP — (J2 _i_ Ki , 3K 2 \ 
^13 - "I- + — + — J 



^14 = 7i (HUT) + IIITI) + IITII) + ITIID) 
E 1 ,= (.J 1 + f + f + §) 

^15 = 7= (HUT) + IITII) - IIITI) - ITIID) 



P — f 7 _L ^2 3A"i , K 2 \ 
^15 - ^-Jl + — — + le ) 



(14) 
(15) 

(16) 
(17) 



Sl ot = -2 



(18) 



^16 = lllll) 

Eie = (Ji + f + f + f) 

We first discuss the ground state (T = 0) entanglement properties. There are 
five distinct eigenvalues: 



ei = E 1 = E 4 = E 9 = E u = E 16 = ( h + f + % + §) 
= E 2 = E 3 = E 6 = E 7 = E 12 = E 13 = — + + ^H 2 ) 

e 3 = E 5 = E 8 = E 15 = (- Jx + f - ^ + §) (19) 



P - P — f 3J 2 , 3Ki , 9K 2 \ 

e 5 = E n = (— 2Ji + f + If + f ) 
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Figure 2: The two resonating valence bond (RVB) states, i[>rvbi( + sign) and 
ipRV B2 ( - sign). A solid line represents a singlet spin configuration. The arrow 
convention is explained in the text 

For simplicity, let us put K\ = K 2 = K. When J 2 < J\ and K < the ground 
state is non-degenerate with eigenvalue e 5 . When J 2 < Ji and K > the 
ground state is three-fold degenerate with eigenvalue e 3 . Thus K = ^ ( J 2 < J\) 
is a QCP. When J\ < J 2 and K < the ground state is non-degenerate 
with eigenvalue e 4 . When J\ < J 2 and X > the ground state is six- fold 
degenerate with eigenvalue e 2 . In this case a QPT occurs at K = ^. With 
K < a QPT occurs at Ji = J 2 when the ground state changes from 
to ^i . In this paper, we focus our attention on this last QPT. The states ipn 
and t/>i describe two resonating valence bond (RVB) states, iJjrvbi and ipRVB2 
respectively. Figure 2 gives a pictorial representation of tPrvbi and ipRVB2- The 
solid lines represent singlets (valence bonds) and the arrow signs follow the phase 
convention that a VB between the sites i and j represents the spin configuration 
^75 [IT (*) 1 U)) ~ II (*) T (j))] 5 if the arrow points away from the site i. 

A measure of entanglement between the spins at sites i and j is given by con- 
currence [HIE]- To calculate this, a knowledge of the reduced density matrix 
is required. This is obtained from the ground state wave function by tracing out 
all the spin degrees of freedom except those of the spins at the sites % and j. Let 
Pij be defined as a matrix in the standard basis {|Tt) > Itl) > lit) > III)} ■ O ne can 
define the spin-reversed density matrix as p = (cr y <S> cr y ) p* (cr y ® a y ), where a y is 
the Pauli matrix. The concurrence C is given by C = max {Ai — A 2 — A 3 — A 4 , 0}, 
where the Aj's are the square roots of the eigenvalues of the matrix pp in descend- 
ing order. C = implies an unentangled state whereas C — 1 corresponds to 
maximum entanglement. The reduced density matrix in the standard basis has 
the structure 



Pij 



(u 



V o 





wi y* 

y u 2 

o o 






v J 



(20) 



and the concurrence CV, has the simple form 
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Cij = 2 max (0, \y\ — ^/wvj (21) 

If the ground state is degenerate, the T = ensemble is described by a density 
matrix which is an equal mixture of contributions from all possible ground states. 
The density matrix is a limiting case of the thermal density matrix as T — > 0. The 
state i^rvbi is the ground state for J 2 < J\ and K < ^r 1 . In this case, the n.n. 
concurrences C\2 = C23 = C34 = C41 = 0.5, i.e., the n.n. spin pairs are entangled 
in equal amounts. The magnitude of the concurrence is independent of Ji, J2 and 
K as long as i/jrvbi remains the ground state. The concurrences C13 and C24 are 
zero, i.e., the spins at the ends of a diagonal are unentangled. At the QCP, J\ = 
J2 = J (K < ^), the concurrences C12, C23, C34, C41 and C^.C^^ are all equal 
to zero. The ground state at this point is doubly degenerate with wave functions 
iPrvbi and ipRVB2- For J 2 > J\ and K < the ground state is given by 
ipRVB2- The n.n. concurrences Cu,C23,Cu and C41 are now zero whereas C13 = 
C24 = 1. The spin configuration described by i\)rvb2 (Fig- 2) can alternatively 
be described as consisting of VBs, i.e., singlets across the diagonals. Since a 
singlet is maximally entangled, C13 = C 2 4 = 1. The entanglement properties of 
a system can further be analyzed in terms of a quantity known as the one-tangle 
Ti which is a measure of the entanglement between a spin and the remainder 
of the system |39[ 1401 141j . It is given as T\ = Adetp^ 1 ' where p« is the single- 
site reduced density matrix. In a translationally invariant system, T\ provides 
a global estimate of the entanglement present whereas the concurrence gives a 
measure of the pairwise entanglement between two individual spins. When T\ = 0, 
there is no entanglement in the ground state, i.e., the state becomes separable. 
The Coffman-Kundu-Wootters (CKW) conjecture [SHI, originally proposed for a 
three-qubit system, can be generalized to yield the inequality 

ri>r 2 = £Cj (22) 

where T\ represents the one-tangle corresponding to the entanglement between 
the ith qubit (spin) and the rest of the system and Cfj is the square of the 
concurrence associated with the pairwise entanglement between the ith and jth 
qubits. The inequality in (22) shows that the pairwise entanglement is not the 
sole entanglement in the system. For the four-spin cluster, t\ has the value 1 
when iPrvbi and/or ipRVB2 are the ground states. The ratio R = ^ quantifies the 
relative contribution of the pairwise entanglement and has values \ and 1 in the 
ground states iPrvbi and ipRVB2 respectively. Roscilde et al. [101 E] have shown 
that the value of R reaches a minimum (not zero) at the QCP of 5* = | XYX 
AFMs in an external magnetic field. In the present case, we have a first order 
QPT. At the transition point Ji=J2=J (K < ^), the ground state is doubly 
degenerate so that the system is in a mixed state. The entanglement measure 7~i, 
defined for pure states, needs to be generalized to the case of mixed states. This 
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is done [39j by considering all possible pure state decompositions of the density 
matrix p. For each of the decompositions, one can determine the average value 
of T\. The minimum of the average over all decompositions is taken to be r{" m 
which replaces T\ in the CKW inequality in Eq. (22). While calculation of r™ n 
is difficult, one can readily see that R [R — at the QPT point either has 

the value zero (t 2 =0, r™™7^0) or is undefined (t 2 =0, r[" m =0) . In the former 
case, the value of R reaches a minimum at the transition point. In both the cases 
R has distinct values on both sides of the transition point. In iPrvbi, two-spin 
entanglements exhaust the one-tangle whereas the opposite is true in the case 
of ipRVB2- A clearer signature of first order QPT is provided by the jumps in 
both the n.n. and diagonal concurrences [TBI [22]. In the present model, the n.n. 
concurrences C12, C23, C34 and C14 are equal to 0.5 in the ground state i/)rvbi and 
zero at the transition point as well as in the ground state ipRVB2- The diagonal 
concurrences C13 and C 24 are equal to 1 in i^rvb2 and zero at the transition point 
as well as in the state iPrvbi- The jumps in the magnitudes of the concurrences 
are associated with the jumps in the density matrix elements, a typical feature 
of first order QPTs HE]. 

We now discuss the finite temperature entanglement properties of the spin 
tetramer. The thermal density matrix, p{T) = j?exp(—j3H) (J3 = ^f), now 
replaces the ground state density matrix with Z denoting the partition function 
of the system. The reduced thermal density matrix pij{T) has the same form as 
in (20) with Cy(T) given by 

Cy(T) = ^max (o, \y(T)\ - yJu(T)v(T)j (23) 
For the four-spin cluster, the thermal density matrix is 

1 16 

P( T ) = y E exp{-/3E k ) |Vfc> M (24) 
z k=i 

where the \ipk)' s and the Ek's are given in equations (3)-(18). The matrix elements 
u, v and y of the reduced thermal density matrix pu(T) are 



y g c 2 c 3 c 



where the eigenvalues e/s (i = 1,2,..., 5) are given in Eq. (19). Due to trans- 
lational invariance, the reduced density matrices for the other n.n. spin pairs 
have the same matrix elements as in the case of pi 2 (T). Figure 3 shows Ci 2 as a 
function of for £ = 0.5 and for r = J {Ki = K 2 = K) = 0.4 (a), 0.2 (6) and 
0.0 (c). Increase in the strength of the four-spin interaction reduces the magni- 
tude of the n.n. concurrence. The value of the concurrence is non-zero provided 
\y\ — y/uv (Eq.(23)) is > 0. One can define a critical temperature T c beyond 
which the entanglement between n.n. spins disappears |S3[42j. One can show 
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0.2 0.4 0.6 0.8 

Figure 3: Concurrence Cu as a function of for ^ = 0.5 and for r = j^(K\ = 
K 2 = K) = 0.4 (a), 0.2 (b) and 0.0 (c). 




0.2 0.4 0.6 0.8 1 

J2/J1 

Figure 4: A plot of -jp-, where T c is the critical entanglement temperature, versus 
J, for r = 0.4 (a), 0.2 (6) and 0.0 (c) 

that in the parameter regime of interest, the thermal entanglement between the 
diagonal spins is zero so that T c can be taken as the critical temperature beyond 
which the entanglement between any two spins is zero. The critical temperature 
T c is obtained from \y\ — y/uv = (Eq. (23)), i.e., as a solution of the equation 

z 3 ~ 3 -¥ - 6z 1+ri+I ? - z 1+2r ^ - 10 = (26) 

where z = e k B T , n = ^ and r 2 = ^. Figure 4 shows a plot of versus ^ for 
r = 0.4 (a), 0.2 (6) and 0.0 (c). For a fixed value of ^, the critical entanglement 
temperature T c decreases as the strength of the four-spin interaction increases. 
T c tends to zero as ^ approaches the QCP ^ = 1. For J 2 > Ji (with K < ^), 
the n.n. concurrences are zero. 

We next calculate the concurrence for pairwise entanglement between the 
spins located at the ends of a diagonal. The matrix elements u, v and y of the 
reduced thermal density matrix pn(T) are given by 
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0.2 0.4 0.6 0.8 1 

k B T/J 2 

Figure 5: Concurrence Ci 3 as a function of for ^ = 0.5 and for r = j^{K\ = 
K 2 = K) = 0.4 (a), 0.2 (b) and 0.0 (c). 



y = | e _/3ei - e _/3e2 + ±e _/3e3 - ^e^ 64 + ^e" 13 " 5 ^ ' 

where the eigenvalues e/s are given in Eq. (19). The reduced density matrix 
P2i{T) has the same elements as in the case of p i3 (T). Figure 5 shows C i3 as 
a function of for ^ = 0.5 and for r = fJyK x = K 2 = K) = 0.4 (a), 0.2 
(6) and 0.0 (c). Again, at a fixed value of the magnitude of C13 decreases 
as the strength of the four-spin interaction increases. The critical entanglement 
temperature T c , beyond which the entanglement between spins located at the 
ends of a diagonal disappears, is also the temperature beyond which the pairwise 
entanglement between any two spins vanishes since in the parameter regime of 
interest the n.n. concurrences are zero at all T. The critical temperature T c is 
obtained solution of the equation 

^i- 2 ? _ 3z 2ri+Z ? - z 3ri ~^ - 5 = (28) 

where z = e k s T , n = j± and r 2 = ^. Figure 6 shows a plot of versus ^ for 
r = 0.4 (a), 0.2 (b) and 0.0 (c). For a fixed value of ^, the critical entanglement 
temperature T c decreases as the strength of the four-spin interaction K increases. 
T c approaches zero as ^ approaches the QCP ^ = 1. The major conclusion one 
arrives at from an examination of Figs. (3)-(6), is that, as in the T = case, 
the two sets of concurrences (i) C 12 , C 23 , C 34 , C41 and (n) Ci 3) ,C 2 4 are mutually 
exclusive. For finite values of the concurrences belonging to the first set, the 
values of the concurrences belonging to the second set are zero and vice versa. 
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0.4 0.6 

Jl/ J2 



Figure 6: A plot of ^|pS where T c is the critical entanglement temperature, versus 
A for r = 0.4 (a), 0.2 (6) and 0.0 (c). 



III. Entanglement Witness 

We now consider the S = ~ polyoxovanadate AFM compound 1/12 and show that 
the magnetic susceptibility x serves as an EW for this compound. The magnetic 
properties of this system are well described by considering only the central square 
of localized V r4+ ions j^E]. These ions form a square plaquette of S — | localized 
spins described by the isotropic Heisenberg AFM Hamiltonian with only n.n. 
interactions (J 2 = K\ = K 2 = in equation (2)). As shown in Ref. [36j, the 
V^12 compound can be treated as a collection of independent S = \ tetramers 
with AFM n.n. interactions of strength ~ 17.6 K. In fact, the theoretical 
expression for the magnetic susceptibility x °f a tetramer gives a good fit (N 
independent tetramers are to be considered in calculating x) to the experimental 
data for V12 (Fig. 2 of Ref. |36j). The susceptibility for a spin tetramer with 
only n.n. interaction of strength J\ is given by 

X 10e~ 3 ^ Jl + 4e~ 2 P Jl + 2e~@ Jl 

(29) 



(gHB) 2 /Ji 1 + e- 2 ^^ + Se-^ + 6e" 2 ^ + Se" 3 ^ 

Following Ref. [28J, the magnetic susceptibility, Xcn along the direction a (a = 
x, y, z) can be written as 

^=H§- < (M * )2 > (3o) 

where M a = J2j denotes the magnetization along a. The expression in 
(30) holds true when the external magnetic field is zero and the Hamiltonian 
is isotropic in spin space. The angular brackets in (30) denote the thermal ex- 
pectation value. The susceptibility Xa can further be written as 
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x* = ^E<W) ( 31 ) 



Due to the isotropy of the Hamiltonian, Xx = Xy = Xz = X an d we can write 

-W [ T + |« (32) 

where N is the total number of interacting spins. The summation of expectation 
values in (32) can be considered as the expectation value of the sum Hs of 
interaction terms describing all-to-all spin couplings. The expectation value of H $ 
has an overall negative contribution to x because of AFM correlations. Hs has 
the nature of a Hamiltonian and the maximum negative expectation value is given 
by the ground state energy of Hs- For separable states, the energy minimum is 
given by the ground state energy of the equivalent classical Hamiltonian [32, 33J. 
For all-to-all spin couplings, the minimum energy separable state is described by 
any spin configuration with total spin vector zero. For N=4 (spin tetramer), the 
classical ground state is given by the Neel state and {Hs) = — §■ For general 
separable states, (H s ) has a lesser negative contribution to x an d one can write 
down the inequality 



{gH B ) 2 2 

x ~ l^rS (33) 

for separable, i.e., unentangled states. Figure 7 shows a plot of ^t^W^ versus 
T (Curve a ) for n independent tetramers, the case of V12. The expression for 
the susceptibility of a single tetramer is given in (29). Curve b represents the x 
versus T curve describing the equality in (33). In plotting the curves, the value of 
jrg is taken as 17.6 K, the experimental estimate for V12. The intersection point 
of the two curves provides an estimate, T c ~ 25.4 K, of the critical entanglement 
temperature below which entanglement is present in V12. The theoretical value 
of the critical temperature, above which the two-spin entanglement disappears is 
obtained from Eq. (26), with n = r 2 = 0, as ~ 15.2 K. Since T c > T C W, 
only multipartite entanglement is present in the thermal state of the tetramer for 
TP <T <T C . 

We now examine whether four-spin entanglement exists in the thermal state of 
the tetramer. This is done by calculating the state preparation fidelity F defined 

as 

F(p) = (^ghz\p(T)\^ ghz ) (34) 

where \i/)ghz) = ^75 (ITITI) + IITIT)) is the four-spin Greenberger-Horne-Zeilinger 
(GHZ) state The sufficient condition for the four-particle (N = 4 ) entan- 



13 



150 
T (K) 



Figure 7: A plot of — — \$ . - ( dimensionless unit ) versus T (Curve a ) for n 

° 1 n(gtiB) /Jl v ' y ' 

independent tetramers, as is the case of V12. Curve b represents the x versus T 
curve describing the equality in (33). The intersection point of the two curves 
represents the critical entanglement temperature T c ~ 25.4 K with ^ ~ 17.6 K. 



glement is given by 

Hp) > \ (35) 

For a tetramer with only n.n. interactions of strength Ji, F(p) is calculated as 

l e -l3Ji _|_ 2 e 2/3Ji 

F(p) = 5e-? J i + 7 + 3e^ + e 2 ^ ^ 

F(p) = |,i.e., > | as T — > indicating the presence of four-spin entanglement in 
the ground state of the tetramer. The critical entanglement temperature, 
beyond which the four-spin entanglement vanishes is obtained from a solution 
of the equation F(p) = \. The value obtained is kB J^ ~ 0.417 from which 
T C W - 7.4 K, assuming j 1 - ~ 17.6 K as in the case of the compound V12. One 
finds that Tjf> is less than the critical temperature ~ 15.2 K. We next 
consider a tetramer with n.n., diagonal and four-spin exchange interactions of 
strength Ji, J 2 and K\ = K 2 = K respectively. Fig. 8 shows a plot of B ^ 
versus ^ for r = ^ = 0.4 (a), 0.2 (b) and 0.0 (c) respectively. For a fixed 
value of ^, the critical temperature for four-spin entanglement decreases as the 
strength of the four-spin interaction increases. 

The tetramer with J\ = J 2 = J and K — 0, alternatively described as the 
tetrahedron, provides an interesting example of the magnetic susceptibility x 
serving as a witness for entanglement other than the entanglement between in- 
dividual spins. The two-spin entanglement vanishes in the thermal state of the 
tetrahedron. The same is true when T = and the system is at the QCP J\ = J 2 - 
Figure 9 shows the EW curves for x (the same inequality bound as in (33) holds 
true) which intersect at a finite temperature ~ 1.9 showing that the thermal 
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Jj/Jl 

Figure 8: Plot of ^J^- versus Jforr = | = 0.4 (a), 0.2 (6) and 0.0 (c). T^is 
the critical entanglement temperature above which the four-spin entanglement is 
zero. 

states are entangled below the critical temperature corresponding to the inter- 
section point. The entanglement is exclusively multipartite in nature. Figure 10 
shows the EW curves for \ in the case of a symmetric trimer described by the 
S — \ Heisenberg AFM Hamiltonian 

Htrimer = J(Si-S2 + S2.S3 + S3. Si) (37) 

In this case, it is well-known [33 ESI that there is no pairwise entanglement both 
at T = and at finite temperatures. In the inequality for \ (Eq.(33)), the factor 
| is replaced by the factor |. In the classical ground state of Hs, the interacting 
spins form angles of ^ with each other. The critical temperature is given by 
^jS- ~ 1.4 . Again, only multipartite entanglement is present in the thermal 
state of the symmetric trimer. 

Another EW, which provides an estimate of critical entanglement temperature 
T c , is based on energy (S21ISH1- The entanglement gap Ge is defined as 

Ge = E sep — Eq (38) 

where E is the ground state energy of the Hamiltonian H describing the system 
and E sep is the minimum energy of the separable states. If Ge is > 0, a finite en- 
ergy range exists over which all states are entangled. For a positive entanglement 
gap Ge > 0, one can define an EW 

Zew = H — E sep I (39) 

where / represents the identity on the full Hilbert space. For any separable state, 
Tr(Z E wPse P )) > 0. If the state is entangled, Tr(Z E wPent)) is < 0. For example, if 
the state belongs to the ground state manifold, Tr(Z E wPo) — E — E sep < 0. Z E w 
thus acts as an EW.The entanglement gap temperature T E is given by U (T E ) = 
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Figure 9: The EW curves for x in the case of a symmetric tetrahedron with 
J\ = <h = J and K = 0. 

E sep , where U(T) (= (H) = — j|§§) is the thermal energy at temperature T. 
For T < T E , the thermal state is entangled and hence T E is a measure of the 
critical entanglement temperature. E sep is given by the ground state energy of 
the corresponding classical spin model [32, 33j. For a square of spins, E sep — —J\, 
as the classical ground state is given by the Neel state. The expression for U(T) 
is obtained from the partition function Z of the square of spins As shown in 
[SHI, in the case of bipartite graphs and lattices, the EW detects only bipartite 
entanglement. Thus Te for the square of spins has an identical magnitude as 
that of Tjp at and above which such entanglement vanishes. In the case of 
non-bipartite graphs and lattices, the EW can detect multipartite entanglement. 
The tetrahedron and the symmetric trimer are examples of non-bipartite graphs. 
E sep in these two cases can readily be calculated as E sep = — 0.5J (tetrahedron) 
and E sep — — | J (symmetric trimer). The entanglement temperature Te has the 
magnitude kB J E ~ 1.9 (tetrahedron) and kB J E ~ 1.4 (symmetric trimer). In both 
the cases, two-spin entanglements are absent and the entanglement present in the 
system for T < Te is multipartite in nature. 

IV. Summary and Discussion 

In this paper, we consider a spin tetramer (S = with n.n., diagonal and four- 
spin AFM exchange interactions of strength J\, J 2 and K x — K 2 = K respectively. 
The significance of the inclusion of three-spin and four-spin interactions in spin 
Hamiltonians of interest has been pointed out earlier (HI EH]. We study the 
ground state and thermal entanglement properties of the tetramer in the various 
limiting cases. At T = 0, QPTs occur as the exchange interaction strengths are 
tuned to certain critical values. We focus on a particular QPT at J\ = J 2 = J 
(K < | J) as the other QPTs exhibit similar features. The QPT point separates 
two RVB ground states, iPrvbi and ipRVB2- The entanglement between two spins 
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Figure 10: The EW curves for x i n the case of the symmetric trimer. 

is determined by calculating the concurrences C12, C23, C34, C41 and C13, C24. The 
n.n. concurrences are non-zero only in ipRVBi and the other two concurrences as- 
sociated with diagonal spins are non-zero only in tpnvB2- The 1— tangle n, a 
measure of the global entanglement, has the value 1 in both iPrvbi and ipRVB2- 
The entanglement measure 7~i, defined for pure states, has to be generalized to 
r{" m at the QPT point where the ground state is doubly degenerate. The entan- 
glement ratio R = has the value zero at the transition point if r{ mn 7^ and 
is undefined otherwise. Away from the transition point, R = 0.5 in the ground 
state ipRVBi and 1.0 in the ground state i/jrvb2- A better evidence of the first 
order QPT is provided by the jumps in both the n.n. and diagonal concurrences 

The study of finite temperature entanglement properties again shows the ex- 
istence of two distinctive parameter regimes. The n.n. concurrences are non-zero 
only when J 2 < J\{K < ^p) and the concurrences associated with diagonal spins 
are non-zero only when J\ < Ji{K < At J\ — J 2 , all the six concurrences 

are zero. The critical entanglement temperature, T c , beyond which entanglement 
between two spins disappears, is computed. The magnitude of T c is highest when 
J 2 = and K — 0. For fixed values of J\ and J 2 , T c decreases as the strength of 
the four-spin interaction increases. A measure of the four-spin entanglement in 
the thermal state of the tetramer is obtained by calculating the fidelity F(p). The 
critical temperature, beyond which the four-spin entanglement disappears 
is calculated and one finds that at fixed values of J\ and J 2 , the magnitude of 
T^> decreases as the strength K of the four-spin interaction increases. 

Molecular or nanomagnets provide examples of spin systems in which the 
dominant exchange interactions are confined to small spin clusters like dimers, 
trimers and tetramers. In several cases, the magnetic properties can be well 
explained by treating the solid to consist of independent spin clusters. We con- 
sider one such compound, V12, which is a collection of spin tetramers with only 
n.n. exchange interactions. Treating the magnetic susceptibility % as an EW, 
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the critical temperature, T c , below which entanglement is present in the sys- 
tem, is estimated from the experimental data on \. The entanglement includes 
both bipartite and multipartite entanglement with T c ~ 25.4 K in the case of 
V12 (-^ ~ 17.6 K). From theoretical calculations, the critical temperature 
beyond which bipartite entanglement vanishes is given by T«~ 15.2 K. bince 
T^ 1 ' < T c , multipartite entanglement in the system persists upto a higher tem- 
perature. The entanglement contents of the thermal states of the tetrahedron 
and the symmetric trimer are shown to be exclusively multipartite in nature. An 
EW based on energy provides evidence of bipartite entanglement in the case of 
a square of spins (relevant for V"12) and multipartite entanglement in the cases 
of the tetrahedron and the symmetric trimer. The EW based on susceptibility 
X can detect both bipartite and multipartite entanglement. The EW based on 
energy detects only bipartite entanglement when the spin system is defined on a 
bipartite graph or lattice. The latter EW can detect multipartite entanglement 
only in the case of a non-bipartite graph or lattice. The critical entanglement 
temperature T c and the entanglement gap temperature T E have identical values 
if the corresponding EWs both detect entanglement. 

After our work was completed, we learnt of a new inequality for susceptibility 
serving as an EW |33I3H|- The inequality can be derived using the sum uncer- 
tainty relation for spin— | operators [33 HHj- When Xx = Xy = Xz = X-> the 
separability criterion for a single cluster of N spins is given by 

The results reported by us, using slightly different arguments, are special cases 
of the general condition (40) for N = 3 and 4. We can generalize our derivation 
in the following manner to obtain (40). We start with the identity 

1 N 

H s = Y,S l S J = -(S 2 -j:sf) (41) 

i<j Z i=l 

where S is the total spin vector. The maximum negative contribution to x is 
obtained for S =0. Thus for separable states with (Sf) = \, Eq.(32) reduces to 
the inequality in (40). There is now a wealth of experimental data on molecular 
magnets and other magnetic systems which are yet to be analyzed in terms of the 
entanglement properties of the systems |35[ 146) . Appropriate finite temperature 
measures of the different types of entanglement need to be developed so that 
contact between theory and experiments can be made. A challenging task ahead 
is to develop suitable EWs which provide signatures of the different types of 
entanglement in the experimental data. 
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